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Formulas for Inverse Osculatory Interpolation 
in the Complex Plane 

Herbert E. Salzer 1 

Improved formulas for inverse osculatory interpolation in the complex plane arc ob- 
tained by inversion of Hermite's formula and the use of appropriate grid point configurations. 
They cover the cases for w = 2(l)7, where n is the number of points required in direct oscu- 
latory interpolation. The formulas provide an improved means for inverse interpolation 
in the complex plane where the first derivative is either tabulated alongside the function 
or is easily obtained. 

Formulas for /i-point inverse osculatory interpolation for finding x=x J rph from /(a) in 
cerms of f k =J(x k ) and f k =f'(x k ), where the Xk^Xo+kh, the k ranging from — [(n— 1)/2] to 
[n/2], are equally spaced at intervals of h, have been obtained previously by the author [1] 2 
through the inversion of the Hermite osculatory interpolation formula [2,3]. Those formulas 
are still applicable to analytic functions when x is replaced by a complex variable z=z -{-Ph, 
and z k = z -\-kh are equally spaced upon any line in the £-plane at interval h, which may be 
any suitably chosen complex number. But for many analytic functions that are given with 
their first derivatives over a Cartesian grid of length h in the complex plane (so that now h 
is a positive number), the accuracy of inverse osculatory interpolation is considerably greater 
when the points z k are chosen as grid points that lie much closer to each other, no longer in a 
straight line beyond the 2-point case. 

Improved formulas for inverse osculatory interpolation, which utilize the Cartesian grid, 
will be particularly useful in connection with such tables for complex arguments as (1) log T(z), 
together with its derivative the psi function, (2) Bessel functions of the first or second kind 
giving Jo(z) and J l (z) = —J' (z), Y (z) and Y 1 (z) = — Yq(z) 1 or linear combinations of them, 
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(3) probability integral I e u "du with its integrand e z [4,8], (4) miscellaneous tables of integrals 

Jo 
of the more elementary functions where the first derivative or the integrand, although not 

tabulated, is easy to calculate, viz., the function (e~ u /u)du [5], and (5) tables of solutions of 

important linear differential equations, together with the first derivative [6]. In all such tables 
and in many others, the user will find the inverse interpolation formulas that are given below 
to be particularly convenient, especially in those cases where grid length h is too large for 
sufficiently accurate complex inverse interpolation, using the formulas in terms of only the 
functional values f t [7]. 

In accordance with the choice of grid points that has been employed consistently in earlier 
papers of the author on complex interpolation, for n ranging from 2 to 7, the w-point formulas 
in this present article are based upon the values of j{z k ) and /'(z*) for z k in the following 
configurations: 
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These formulas enable one to find z=z +Ph by expressing P in terms of j{z^ J rPh) = /, 
J(z k )=J k , &i\df'(z k ) =f' k , for ^ = s +^, where k ranges over the points of the selected w-point 
configuration, the choice of n depending upon the number of points required in direct osculatory 
interpolation. Although the direct interpolation formula for n=6 and n=7 is of the 11th- 
and 13th-degree accuracy, respectively, (recalling that the n-point direct osculatory interpola- 
tion formula is of (2n— l)th-degree accuracy), the inversion formula for P given below does 
not go beyond the lOth-degree terms. In fact, most practical problems will require only the 
first few terms, and it is very rarely that one would need this inversion series up to the 10th- 
degree terms. Thus in this paper we do not use the coefficient of P n in the 6-point Hermite 
formula, nor the coefficient of P 11 , P 12 , and P 13 in the 7-point Hermite formula. 

For every n we define r= (f—ftyj/hfo, and corresponding to each of the above ?i-point con- 
figurations, where n=2(l)7, quantities s, t, u, v, w, x, y, z, and z are defined as follows: 

Two-point 

«=(-3{/o-yi}-A{2/S+/I}W5, 

«=(2{/ -/i}+A{/s+y{})/*/o, 

u=v=w=x=y = z=lz = 0. 

Three-point 

•=*(&yo+(i-4ty 1 -(i+4i)/ i +A{-(4-4t)y5+*yi-/;})/A/i, 

t = U-a2 + l2i)f +(7 + 5i)J l + (5+7i)f i +h{-Sif (> -(2+i)f l + (2-i)f i })/hf , 
u=UWo-(7-2i)fi-(7+2i)f i +h{(4+4i)f +(2-i)f[-(l-2i)f t })/hf , 

»=i(-(4-4»y,+ (l-3Cl/i+(3-*Vi+A{-2fi+^--*y{})/*/5, 

w=x=y=z=~z = 0. 

Four-point 

s=l(l^/o-(8+6i)/ 1 + (8-6t)/ i -7ty i+i +A{-(6-6t)/o+2/;-2 l y;-(l-i)/;+ i })/¥;, 

< = i(-(37 + 37t)/o+(32-12t)/ 1 -(12-32z)/ i +(17+17t)/ 1 +i+A{-17ty;-(4-6i)/;+(4 + 6f)/; 

— 5ifi +i })lhJ'o, 
«=K7Q/o-(15-50tl/ 1 -(15 + 50i)/ i -40/ 1+i +A{(14 + 14f)/;-(4+lH)/i-(ll+40/; 

+(6+6t)y; +i })/¥o, 

»=!(- (38-38»U- (25 + 370./, + (37+25i)f t + (26-26i)f l+i +h{-Uf + (10+3*1/1 
+ (.l0-3i)ff-8fi +t })/hf 0> 

w=4(-23*y +(21+2«U-(21--2»V < +18tyi + .+A{(4-4tVi-(4-3iOyi-(3-4«)/i 

+ (3 -3i)y ;+,.}«, 
x=K(3+3oy -(3-3i)y 1 +(3-3f)y i -(3+3i)y 1+i +A{iy;-tyi-ty;+i/; +i })/¥o, 

Five-point 

s=(^{-(1875-6250t)y - (4000 + 6000t}/i + (103+254tVa + (2272+496^ 
+ (3500- 1000i)f l+i ] + A{ - (200-150^+200^- (3+4i$+ (32-24ty« 

-(5o+50i)y; +t .})Ay; 
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t=(ih{ - (1525+3950*)/ + (5600+ 1200*)/- (223 + 114*)/,- (1952-464*)/, 

- (1900-2400*)/i+,} +M (65-230*)/;- (240-240*)/;+ (9+2*)/ 2 - (32-64*)/$ 

+(i2o+2o*y 1+4 })A/; 

M=(To 1 oo{(59750+37125i)/o-(96000-53000i)/ I +(5284-1313i)/ 2 +(29216-26312i)/ i 

+ (1750-62500*)/+,} + T ^{ (625+1625*)/;+ (100-3400*)/;- (89 -48*")/ 2 - (24+732*)/, 

-(1125-575*)/ 1+< }W; 
»= Gnnn> { - (72625— 5250*)/ + (48000-118000*)/ - (4103-5746*')/ 2 - (14272-43504*)/, 

+ (43000+63500*)/+,} +4o{-(2800 + 1825*)/;+(5200 + 6200*)/;+(96-247*)/ 2 

+ (896+1528*)/;+(1450-2800*)/i +< })/A/; 

w=( T ^ T f{(3850-3075*)/o+(2450 + 9200*)/ 1 -(42 + 581*)/2-(658 + 3344*)/ j -(5600 + 2200*)/i+«} 
+A{(215-5*)/;-(600+110*)/i+(7+21*)/ 2 -(114+72*)/;+(30+250*)/ 1+< })/% 

*=(Fb{- (3500- 10125*')/ -(17250 + 14750*)/i + (1022 + 1121*)/ 2 + (4978 + 6254*)/, 

+ (14750-2750*)/i + ,}+ T ^ r {-(375-250*')/;+(1300-650*)/ l -(49+32*)/ 2 +(316+38*)/; 

- (400+500*)/;+,}) A& 

y=( T TrW{ — (1000 + 5375*)/o+(12500 + 2000*')/ 1 - (898 + 189*)/2- (3852 + 1936*)/, 
-(6750- 5500*)/+,} + T ^{ (50- 100*)/;- (200-350*)/;+ (18-*)/ 2 - (82-24*)/; 
+ (175+75*)/;+,}) A/ 0J 

*=(tAt{ (375 + 500*)/ - (1500-500*)/ + (117-44*)/ 2 + (508+44*)/,+ (500-1000*)/+,} 
+4o (25*/;- 100*/;- (4 - 3*)/ 2 + (16 - 120/;- 50/;+,} )/hf , 

2=0. 

Six-point 

s=(ttJW{16500*/ - (14336+ 1152*)/ + (289-98*)/ 2 + (14336-1152*)/,- 14000*/, +< 

- (289 + 98*)/ 2 ,} + T ^{ - (400-400*)/;+ (192-256*)/;- (4-3*)/ 2 + (256-192*)/$ 

- (200-200*)/;+,- (3-4*)/ 2 ',})//7/;, 

<=(innnr{ — (86500+86500*)/o+(91392 — 102656*)/i— (1219— 2983*')/ 2 — (102656-91392*)/* 
+ (96000 + 96000*')/ i +,+ (2983-1219*)/ 2 ,} + 2 7t{-3100*/; + (896+3072*)/;+(4-53*)/2 

- (896-3072*)/;-2800*/;+ t - (4+53*)/ 2 ,})/A/;, 

«=(oto{95150/ + (18688+ 124416*)/- (1463+2934*)/ 2 + (18688-124416*')/—129600/ 1+i 

-(1463-2934*)/ 2 ,}+ ¥ \{ (1500 + 1500*)/;- (2816 + 1152*)/;+ (38+39*)/ 2 - (1152+2816*)/; 
+ (1920 + 1920*)/;+,+ (39 + 38*)/2,})//\/;, 
«= (Wirol — (169500- 169500*)/ - (301824 + 192768*')/! + (9957 +2601*)/ 2 + (192768 + 301824*)/, 
+ (271200-271200*)/ 1 +,-(2601+9957*)/ 2i }+ T | T y{-4970/;+(7808-4224*)/;-(172-9*')/; 
+ (7808+4224*)/;-8160/;+,-(172 + 9*)/;,})///;, 
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^=(8o 1 oo{- 104 7750iy +(1662976-466368O/i- (44774-307430/- (1662976+4663680/* 
+ 1919000*/ 1 +*+(44774+307430/ 2 *} + sio {(14650- 146500/J- (8672-424960/1 
+ (564-7230/2- (42496-86720/*+ (29300-293000/1+*+ (723-5640/^} Wo, 

*=(T^^{(919000 + 919000i)/ -(1038592--20802560/i + (8369-79033i)/2 

+ (2080256- 10385920/,- (1890000 + 1890000i)/i + i—(79033-8369i)/2f}+TArr{ 2 4850zyi 

- (36096+462720/1 + (271+ 13280/2 + (36096-46272O/J + 58600i/; +i 

-(271-13280/ 2 *})/¥o, 
y=(sio{ -28450/o-(19568 + 509760/i + (1393 + 14490/- (19568- 509760/*+ 64800/+* 

+ (1393- 14490/2*} +t4o { - (750 + 7500/;+ (2752 + 1440/1- (61 +33z)/£+ (144 +27520/ 

-(2040 + 20400/1+ — (33 + 6l0/a*})/A/n, 
3-( T ^oo{ (12000- 120000/o + (32224 + 127680/- (1357-990/- (12768 + 322240/z 

-(30000-300000/+/- (99- 13570/ 2 z}+ T |o{310/;- (608-6240/1+ (22-90/2 

-(608 + 6240/; + 960/l +( -+ (22 + 9i)f' 2t })/hf , 
5= (8oVo{15750i/ - (30688- 147840/ + (812-11590/ + (30688 + 147840/z-43000z/ 1+i 

-(812 + 11590/ 2 a + s^o{ -(200-2000/;- (64 + 8480/1- (7-240/ 2 + (848 + 640/; 

-(700-7000/;+*- (24-70/2 J)/A/o. 

Seven-point 

*=(rinnr{-(1960-217800/o- (38400- II2OOO/1- (1119 + 13420/2+ (15744 + 145920/* 

+ (28000-460000/+*-(2432-640/ 2 +z+(167-2940/ 2 *}+t^o{-(480-4400/o 
-800^1+(24 + 7O/ 2 + (384 + 112O/;-(1400 + 2000/l +i +(32 + 16O/ 2 +*-(3+4O/ 2 *})/^ 
t=(-2Woo{- (H3316 + 1470120/o+ (193280-4070400/ + (21313 + 16590/ 2 - (242688 + 307840/* 
+ (120000+6000000/i+*+ (20736-209280/ 2 +*+ (675+41050/ 2 *} +-s-U{(344-40920/i 
+ (7040+60800/1- (268- 1510/2- (4448- 15360/*+ (13600-92000/1+*- (416-1280/ 2 +* 
+ (60+150/ 2 *})/¥o, 
u= (Woo { (825894 + 140208O/o+ (898176 + 22840320/- (94875-832500/ 2 
+ (1152000-6720000/*- (2760000 + 19920000/+*- (2688- 1728000/ 2 +* 

- (18507 + 16290O/ 2 *} +4o{(10044 + 129080/o- (51392-66560/1+ (470-17650/ 2 
+ (14880-217600/;- (21600-904000/;+* + (1248-22400/ 2 +*- (333- 1580/2*})A/^ 

*= (w^o { - (2094204 - 13613720/o- (7780992 + 32285440/ + (24441 -4663870/ 

- (1650816-42565120/*+ (11844000- 14600000/+*- (433920+4663040/ 2 +* 

+ (91491+335l0/ 2 *} + 8w{-(44834 + 74880/ +(103744-1473920/; + (3484 + 58370/2 
+ (4704 + 918720/*- (160800+2896000/;+*+(2400 + 91840/ 2 +*+(588-12450/ 2 *})/% 
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w= (soVo {(662260-25734300/o+ (9236480-39006400/ + (404714 +4376770/ 2 
-(1662464+51515520/i— (9358000— 11071000f)/i+i+(810752+46496f)y2+< 
-(93742-704490/ 2 i}+^{(38730-25190^+(43840+2112800/;- (8494+20170/2 
-(71904+83672iy«+(387500+130100^+«- (10352+^^ 

»=(T^^{(1955376 + 36844320/o-(6505920-15994560?:)/ 1 -(1095693 + 27999i)/2 
■+(7805568 + 5022624i)/ i -(1074000 + 25434000i)/i +f -(1139136-969888i)/ 2+l - 
+ (53805-209505zy 2i }+ T ^ nr {-(17684-67962f)/i-(315360+197120i)/ 1 
+ (13623-86860/2 + (186528 + 371040/J- (667600-2930000/1+*+ (22016 -46880/ 2 +* 
-(2615+20100/ 2 *} Wo, 

V= (Woo {-(H49156 + 4146420/o-(2086704 + 50561280/i + (264015-2577300/ 2 
-(2906640-3424800/i + (5796000 + 58230000/i+<+(52992-5037(H)0/2+f 
+ (29493 + 667800/2i}+^w{- (9306 + 177920/o+ (1141 68-278240/1- (1095-53650/2 
-(53460-273200/*+(98100-2179000/l + *-(4332-63200/ 2 +* + (1122-120/2f})A/o, 

0=( T ^ Tnr {(985O56-4OO8O8O/o+(4675648+2O23936O/i+ (221+3634030/ 2 
+ (1877904-19585280/— (7782000 + 3940000/i +i + (306480 + 4002560/ 2+l 
- (63309 + 342590/ 2i } +ttoo { (16076 + 59820/J- (56336 -932480/i- (3071 +44780/ 2 
+ (23024-505680/*+(69600+2194000/i + i-(1200+75760/ 2 +i--(877-6900/2<})A/o J 

i=(W^{-(74180-1509900/o-(802080-2974400/i-(45912+448910/ 2 -(35088-4520160/i 
+ (1034000-8390000/ w -(89296 + 146080/ 2+ t+(12556-19470/2*}+^ir{~(2490-9700/o 
-(4240+176800/i+(952+16lO/2+(2232+89760/*-(34300+205000/i + z+(1096+8480/2^ i 
+ (51-1920/ 2 a)/¥o. 

For every n, P is given by the following formula: 
P=r-r 2 s+rH2s 2 -t)+r\-5s* + 5st-u)+r 5 (Us±-2lsH + 'M 2 + 6su 

-28s 2 u+7tu + 7sv— w)+r 7 (132s 6 — 330sH+l80sH 2 +120s 3 u— 12/ 3 -72.^— 36s 2 v+4u 2 +8tv 

+ 8sw-x)+r 8 (-±29s 7 + 1287s 5 *- 990s¥-495s%+495s 2 ^+165sf 3 + W5s 3 v-45t 2 u—4:5su 2 

-90stv-A5s 2 w + 9uv + 9tw+9sx—y)+r 9 (U30s s — 5005s 6 £+5005s¥+2002s 5 ^— 1430s 2 * 3 

-2860s 3 ta-715s 4 i?+55* 4 +660$Z 2 M+330s% 2 +6 

-1108^-55s 2 ;r + 5?; 2 + 10^/^+10^+10^— O+^ 10 (— 4862s 9 + 19448s 7 *— 24024s 5 * 2 — 8008s 6 ?/ 

+ 10010s¥+15015s 4 *u+3003s 5 ^-1001s* 4 -6006s%^ 

+ 286* 3 ^ + 858s*^ 2 + 858s 2 *w + 858s* 2 z; + 858s 2 ^ + 286^ 

-132suw—lS2stx— 66s 2 y+nvw+llux + llty+llsz-'z) + . . . . 
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